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Abstract 

In this paper we analyze the relativistic quantum motion of charged spin-0 and spin-i particles 
in the presence of a uniform magnetic field and scalar potentials in the cosmic string spacetime. In 
' order to develop this analysis, we assume that the magnetic field is parallel to the string and the scalar 

potentials present a cylindrical symmetry with their center on the string. Two distinct configurations 
for the scalar potential, S(r), are considered: (i) the potential proportional to the inverse of the 
. polar distance, i.e., S oc 1/r, and (ii) the potential proportional to this distance, i.e., S oc r. The 

J> " energy spectra are explicitly computed for different physical situations and their dependences on the 

^vO . magnetic field strength and scalar coupling constants are presented. 

OC 1 PACS numbers: 03.65.Pm, 03.65.Ge 

cn 

! 1. Introduction 

O 
_ _ . 

The analysis of a quantum mechanical system under the influence of the gravitational field has attracted 
attention in particle physics since many years ago. In this direction, the hydrogen atom in curved 
spacetime has been investigated by IT]|30j|3l. It was shown in [30] that the shifts in the energy spectrum 
caused by local curvature is different from the usual gravitational Doppler shift. However, this shift 
is appreciable only in the region of strong gravitational field. Recently, the analysis of the influence 
of the spacetime produced by topological defects on the energy spectrum of the hydrogen atom, has 
been considered under both non-relativistic JH and relativistic [5] points of view. In these papers, the 
hydrogen atom is placed in the spacetime produced by an idealized linear cosmic string and a point-like 
global monopole. Different from the previous analysis, due to the specific geometry associated with the 
spacetime produced by these idealized topological defects, in these more recent investigations the energy 
spectra associated with the hydrogen atom can be exactly obtained. 

According to the modern concepts of theoretical physics, different types of topological defects 
may be produced by the vacuum phase transition in the early Universe (UI3. These include domain 
walls, cosmic strings and monopoles. Among them, cosmic strings and monopoles seem to be the best 
candidates to be observed. The former are linear defects, and the spacetime produced by an idealized 
cosmic string is locally flat, however globally conical, with a planar angle deficit determined by the string 
tension. Due to this conical structure, a charged particle placed at rest in the cosmic string spacetime 
becomes subjected to a repulsive electrostatic self-interaction ||8] HI proportional to the inverse of the 
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polar distance from the string. Also it was shown in |[TU1 that linear electric and magnetic sources in the 
spacetime of a cosmic string become subject to induced self-interactionsQ 

The non relativistic quantum mechanical analysis of the motion of a charged particle in the cosmic 
string spacetime taking into account the electrostatic self-interaction and in the presence of a uniform 
magnetic field parallel to the string has been developed in IfTD . There it was shown that the presence 
of the defect reduces the degeneracy degree of the Landau levels. In addition, the inclusion of the 
electrostatic self-interaction in this analysis breaks the residual degeneracy. In both situations, exact 
wave functions and energy eigenvalues have been found for specific magnitudes of the magnetic field. 
Moreover, in 11121 the analysis of the non-relativistic motion of two charged particles on a cone and in 
the presence of a static magnetic field has been developed, too. 

The well-known procedure to introduce the coupling between a charged particle and electromagnetic 
fields in the Dirac and Klein-Gordon equations, is through the minimal coupling. By this procedure, the 
four-vector differential operator - id^ is modified in order to include the electromagnetic four-vector, 
Ap, as Pn — > p^ - eA M . Some years ago, Dosch, Jansen and Miiller in lfl3l . pointed out that the minimal 
coupling is not the only way to couple a potential to the Dirac equation. There, it was suggested that a 
non-electromagnetic potential can be taken into account by making a modification in the mass term as 
M — > M + S(r, t), being S(f, t) the scalar potential. This new formalism has been used by Soff, Miiller, 
Rafelski and Greiner in lTT4l . to analyze the Dirac equation in the presence of a Coulomb potential and 
a static scalar potential proportional to the inverse of the radial distance. More recently, the electrostatic 
self-interaction on a charged particle in the cosmic string spacetime has been considered in [15 ] and 
lfl6l . as a scalar potential to analyze its relativistic quantum motion. So, if one wants to investigate the 
relativistic quantum motion of a charged particle in the presence of electromagnetic and scalar potentials, 
both procedures, the minimal coupling and a modification in the mass term, should be taken into account. 

In this paper we shall analyze the relativistic quantum motion of bosonic and fermionic massive 
charged particles in the spacetime surrounding an idealized cosmic string in the presence of a uniform 
magnetic field parallel to the string and also in the presence of scalar potentials. Two specific 
configurations to the scalar potential will be considered: (i) the potential proportional to the inverse 
of the polar distance to the string, S oc 1 fr, and (ii) the potential linear to the polar distance to the string, 
S oc r. We shall see that, for particular values of the magnetic field, we are able to obtain exact solutions 
to the Klein-Gordon and Dirac equations, providing the energy spectra associated with bound states. 

Spherically symmetric linear scalar potential has been introduced into the Dirac equation in ifTTl as a 
simple model to confine relativistic quark states. In this paper some analytical results could be compared 
with those provide by the bag model. In addition to the linear scalar potential a scalar Coulomb-like 
potential has also been included in [18]. In the latter, separete analysis of the energy spectrum are 
examined for the case of absence of the confining potential, followed by the complete analysis, i.e., 
taking into account the linear potential. In these two papers, due to the presence of linear potential, only 
approximated solutions for the energy spectra have been obtained. The spectra found are similar to the 
harmonic oscilator. Moreover, exact solutions to the Dirac equation in presence of an uniform electric 
field and a linear scalar potential have been found in |[T9l . In the latter the zero-component of the vector 
potential and the linear scalar potential are both proportional to the z-coordinate. 

This paper is organized as follows: In section [2] we introduce the physical scenario, namely the 
explicit expressions for the metric tensor, for the four-vector potential and for the scalar potentials. We 
consider these background fields in the analyses of relativistic quantum motion of charged bosonic and 
fermionic particles. The first part of this paper deals with the bosonic case. In section [33 the modified 
Klein-Gordon equation is presented. Considering stationary states, v P(r, i) = e~ lEtx ¥s(f), and admitting 
for an ansatz, that obeys a cylindrical symmetry, we provide, in subsection I3..1I the most general 

differential equation obeyed by the radial function. The obtainment of solutions for this differential 
equation is considered in section [4] specifying separately in each subsection, different values attributed 

'These induced self-interactions are consequence of a distortion on the electric and magnetic fields caused by the planar 
angle deficit produced by the defect. 
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to the parameters associated with the magnetic field and scalar coupling constants. In the second part, we 
repeat this analysis for the fermionic case. In section [53 we present the modified Dirac equation, and our 
choice for the curved Dirac matrices. The following section deals with the radial differential equations 
resulting after the consideration of an ansatz suggested by the symmetry of the system. Solutions to these 
equations and physical informations are obtained in section |T] In section [T] we summarize all results. 
Finally in the appendix [A] we present some properties associated with the most relevant differential 
equation presented in this paper, that is the biconfluent Heun (BCH) differential equation. Here we use 
natural units, % - c - G = 1. 



2. The System 

In this section we present the physical environment under which we shall consider relativistic charged 
particles. We first introduce the spacetime background of an idealized cosmic string where the particle 
will propagate, followed by vector potential associated with a uniform magnetic field parallel to the string 
and the most general scalar potential which allow us to obtain analytical solutions of the corresponding 
Klein-Gordon and Dirac equations. 

By using cylindrical coordinates with the cosmic string along the z-axis, the corresponding metric 
tensor is defined by the line element below: 

ds 2 = dt 2 - dr 2 - a 2 r 2 d<p 2 - dz 2 , (1) 

where the polar coordinate, r > 0, the coordinates (?, z) e (-°°, oo) and the angular variable cf> e [0, 2n\. 
The parameter a, smaller than unity, is given in terms of the linear mass density of the string, fi, by 
a = 1 - 4ju. 

In the cosmic string spacetime, the vector potential associated with an uniform magnetic field parallel 
to the string, B = V x A = Bq k, in Coulomb gauge can be expressed by, 

A = (0, A , 0), with A = - l -aB Q r 2 . (2) 

Also we shall consider in this analysis the presence of cylindrically symmetric scalar potentials which 
allows us to provide analytical solutions for the field equations. This scalar potential readsH 

S( r ) = H£ +T]L r. (3) 
r 

In the next sections we shall develop the relativistic quantum mechanical analysis associated with 
spin-0 and spin- 1/2 charged particles under the influence of the above field configurations. We shall 
start with the bosonic case, followed by the fermionic one. 



3. Bosonic Fields: Klein-Gordon Equation 

The relativistic quantum motion of a charged bosonic particle in curved space, and in the presence of 
electromagnetic and scalar potentials, is described by the modified Klein-Gordon equation below, 

[£> 2 - £R + (Af + S (r)) 2 ] OO) = , (4) 

where the differential operator above is given by 

n 2 = ^=D,(^- gg ^D v ) , (5) 

with = dft + ieA M , g = det^y) and S(r) is the scalar potential. Since we are considering a minimal 
coupling, the previously introduced non-minimal coupling £ needs to be discarded. 

2 The relativistic quantum analysis of the motion of a charged spin-0 particle in the spacetime produced by an idealized 
cosmic string and in the presence of magnetic field produced by a magnetic monopole and scalar potentials were developed in 
(20). 
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3..1 Bosonic Radial Equation 



Here in this subsection we shall construct the radial differential equation associated with the Klein- 
Gordon equation (@]), considering the most general physical situation involving the presence of the 
magnetic field and the scalar potentials. 

Due to the cylindrical symmetry of such system, let us adopt for the wave function the ansatz below: 

<B(jt) = e iikz+m *- Et) R{r) . (6) 

Substituting © into ©, we find that the radial function, R(r), must obey the differential equation 



d 2 ld_ 

dr 2 r dr 



R(r) + 



E 2 -k 2 + 2mJ- - M 2 to 2 r 2 
a 



nr 



a 2 r 2 



-(M + S(r)T 



R(r) = , 



where 



co - 



eB Q 
2M ' 



(7) 



(8) 



is the so called cyclotron frequency of the particle. Considering the scalar potential ©, the equation © 
is expressed as a biconfluent Heun differential equation ETTl . 

In order to solve the above differential equation, it is convenient to define a new radial function, u(r), 
as shown below: 



u(r) = \f?R(r) . 

Taking the scalar potential (f3]) and defining a new dimensionless variable 

x = VXr with A - ■Jrfa + M 2 a> 2 , 
the radial biconfluent Heun equation fTJ takes on the following Schrodinger form 



where 



d u(x) 
dx 2 

7m 



eic,m Jm SlX 2 
— r- + — T "t" + — X 



u(x) = 



E 2 - k 2 - M 2 + 2 (moj- - 7] C ri L ^ , 



m 



2M 2M 
—t]l and 8 C —7]c ■ 

Va Va 



(9) 
(10) 

(11) 

(12) 
(13) 
(14) 



After we have constructed the complete radial differential equation, the next steps are to obtain the 
complete wave functions set and the corresponding energy spectra. 



4. Analytical Solutions for the Klein-Gordon Equation 

In this section, we shall develop the analysis of the radial differential equation obtained in the previous 
section considering different physical situations. For completeness only, we start with the relativistic 
version of the Landau levels for bosonic charged particles. For this case the solutions for the radial 
differential equation are presented in terms of the confluent hypergeometric function. The newest and 
most relevant analysis about the relativistic quantum motion contained in this paper are exhibited in the 
following subsections. There we shall see that the solutions for the radial equation are given in terms of 
biconfluent Heun functions. 
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4..1 Landau Levels 



The simplest particular solution to (fTTT) can be found by admitting rjc = Vl = but keeping co + 0. As 
we know, only bound states exist in such condition. The differential equation (TTTT> can be rewritten as 
follows: 



d u 

dx 2 



j r? 1 1 m 2 

Az 2m J ~ , 
+ — + ^ x 



0. (15) 



It is helpful to analyze the asymptotic behavior of solutions to the equation above, for x — > and 
x — > oo. This analysis indicates that it is useful to write u(x) in the form 

u(x) = Je~^F(x), (16) 

where p = ± + ^,andF(;t) is an auxiliary function. In fact, substituting ([TBI into the differential equation 
we get, 

d 2 F / . . lS dF (E 2 -k 2 -M 2 2m 



x- 



i ~ lX dF (E 2 -k 2 -M 2 2m - \ 
+ (2/8 - 2x 2 ) — + + 2B - 1 )xF = . (17) 



c?x 2 

Defining a dimensionless variable y - x 2 - Mtor 2 , the above equation yields 



d 2 F (y) dF(y) 
y , 7 +(c-y)— aF(y) = 0, (18) 

with the constants being given by 

\m\ l(\m\-m \ E 2 - k 2 - M 2 

c=— + l,a=-\— + 1 . (19) 

a 2\ a J AMoj 

The solution to the differential equation (fT8l ) is the well-known confluent hypergeometric function: 

/1/ N-m \ £ 2 - ^ 2 - M 2 |m| \ 



The asymptotic behavior of a confluent hypergeometric function for large values of its argument is E51 : 

iFKa, b, z) * ^* [i + 0(12)-!)] . (21) 

T(a) L J 

So, due to the divergent behavior of the function F(y) for large values of its argument, bound states 
solutions can only be obtained by imposing that this function becomes a polynomial of degree n. In this 
case, the radial solution presents an acceptable behavior at infinity. This condition is obtained by 

l(\m\-m \ E 2 -k 2 -M 2 

+ 1 = -n, n = 0,1,2,3... (22) 



a J AMoj 

This equation provides the quantization condition on the energy spectrum of the particle: 



E k , m n = ±\k 2 + M 2 + 2Moj 



\m\ — m 
2n+l + — 

a 



(23) 



The presence of the factor a visibly modifies the degenerate spectrum of the particle. Each set of wave 
functions with the same value of n is called a Landau level. In the case of a = 1, the non-relativistic 
analyses of Landau levels is most developed by using Cartesian coordinates, and the spectrum can be find 
in many textbooks. However, by using cylindrical coordinates system, the non-relativistic analysis has 
been developed in [22l as a separate problem and also in |[23l ; as to the relativistic approach this problem 
has been considered in [24J. In [22J the approach adopted to find the energy spectrum is similar to ours; 
there it is shown that the energy depends on the angular quantum number in the same combination as 
presented in the equation above! 



3 Also in 1111 the energy levels depend on the angular quantum number in the same manner. 
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4..2 Vanishing magnetic field 



Let us now admit that oj - 0, but take into account the effects of the scalar potential, i.e., keeping 
non-vanishing 77c and t]l- In this case we can simply rewrite equation (fTTb as followsQ 



d u(x) 
dx 2 



h 7m S C S L X 2 
H — H 1 X 



\rj L \ 



\riL\ 



u{x) = 0, 



(24) 



with 



5z 



E 2 -k 2 -M 2 - 2j]cJ]l , 

2M - 2M 
-_t]l and 6 C = i=Vc 



(25) 



In order to obtain the energy spectrum to equation (l24l ). it is convenient to analyze its asymptotic behavior 
for x — > and x — > 00. Through this analysis it is possible to express the function u(x) in terms of the 
unknown function F(x) as shown below: 



i(x) = xPe 



F(x), 



with 



m- 



Now, substituting the ansatz (1261 ) into (1241 . the resulting equation reads: 



xF"(x) + 1 2/3 + — jc - 2x 2 F'{x) + 



\Vl\ \\Vl\ 4rf L 



F(x) = . 



(26) 



(27) 



(28) 



This equation is the biconfluent Heun's differential equation Ell , whose solution is the so called 
biconfluent Heun (BCH) function, Hg- 



<5l e* 



F(x) = H B 2/3-1, 7^7 + -^-, 25c, 

toil toil W L 



We shall use the Frobenius method to obtain solutions for 

00 

F{x) = ^ c n x" 

n=0 



(29) 



(30) 



Substituting the series above into d28l . the coefficients are defined by a three-term recurrence relation, 

1 



Cn+2 



(n + 2)(2/3 + n + 1) 
and by the relation below, 



5l -\ 
-■ — -{n + 1) - B \c n+ \ + (2n - C)c n 



tod 



/ gi + nc_ 
VtoZfW 1 /? 



c , 



(31) 



(32) 



where£=/3^ + 5 c andC = ^-2/3-l + {|. 



4 The case with a> = rj L = has been analyzed in 1151 . 
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A special kind of exact solutions representing bound states, can be obtained by searching for 
polynomial expressions to F{x). Being n the order of such polynomial, two different conditions must 
be simultaneously satisfied: 



C = In 



(33) 



being n a positive integer number, and 



= . 



(34) 



Admitting that cq = 1 and using the recurrence relation (DTI ) we present below the first coefficients for 
the Frobenius series: 



1 



C2 



C3 



2(2/3 + 1) 
M 



2M 2 



1L . VC 



\\t}l\ 
1 



P 



rV^- toe + 7^0 + 1) -c 



\m 



\Vl\ 



2M 1 



\Vl\ \\Vl\ fi 



(35) 



The condition C -2n provides the following energy spectrum for the particle: 



Ekr 



k z + 2r] C r} L + 2\Tj L \\n+p+- 



(36) 



This expression holds for both signs, and the quantity inside the brackets is clearly positive if tjl and r\c 
have the same sign, so the system allows positive and negative solutions for the energy spectrum. At 
first glance, this equation indicates that the energy spectrum does not depend on the mass of the particle. 
However, we would like to mention that that the condition c n+ i = 0, provides an algebraic expression 
involving specific values of the scalar potential coupling constant, the mass of the particle and angular 
momentum quantum number. Let us assume that the parameter tjl can be adjusted to fit this equation. In 
this case, this parameter depends on the quantum numbers n and m, which we define now by t]l„ m ■ 

In what follows, we shall consider for the function F(x) the simplest case only, a polynomial of first 
order. So, after a few steps we can write: 



F ljn (r) = \ + M[l + U£\ K 



(37) 



Also for this specific case the energy reads, 



Ek,m,\ ~ — 



k' + 2rj c r lL ^ + 2\r lLl J\p+- 



(38) 



Moreover, the condition (T34l) allows us to express the parameter 77^, m in terms of the mass of the particle 
as shown below: 



M~ 

nu, m = -r[nc + W +Dric+P(P+ 1)] • 



(39) 



So the corresponding energy is completely defined, and we can say that the energy spectrum of the 
particle, in fact, does depend on its mass, as should be. Moreover, the energy spectrum also depends on 
the parameter a through the definition of 
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4.3 Linear Confinement 



In this subsection we shall analyze the situation where the the Coulomb-type scalar interaction is absent, 
but keeping the linear confining potential and a magnetic field. For this system, the function u{x) satisfies 
the differential equation below, obtained directly from (fTTT >: 



d u(x) 
dx 2 



u(x) = , 



where 



799 m 1 m 
ek m = E - k - M +2Mu>— and f m = - 

a A a 1 



(40) 



(41) 



Again, analyzing the asymptotic behavior of the differential equation for large and small values of the 
variable we verify that u{x) can be expressed by means of an unknown function G(x) as follows: 



where 



u{x) = A~H-*-x)g(.*) , 



2 a 



(42) 



(43) 



Substituting (l42l into (1401) . we can see that G(x) satisfies the BCH differential equation with one of its 
parameters being zero: 



x G"{x) + {if} +^x- 2x 2 ) G'{x) + 
Consequently, we can write 

G{x) = H B 
Again, using the Frobenius expansion, 

oo 

G{x) = Y J d n x n , 

into (1441) . we find that the coefficients of this series obey the following recurrence relations: 



€ Or 

- + — - 2/3 - 1 I x 
A 4A 2 1 



G(x) = 



2—, — , - + — 0, -x 

a A A 4A 2 



(44) 



(45) 



(46) 



d n +2 - 



1 



(n + 2)(2/3 + n + 1) 



-—(p + n+ l) d n +\ + (2n - C)d n 



and d\ = -—do 
2A 



(47) 



where C = ^- 2/3-1 + jh, and we have set do equal to unity. Below we provide two additional 
coefficients: 



di = — 

d 3 



1 



2/3+ 1 
1 



Mco 3 



VL 20 + 2) 



608+ l)dZ>VM^ 2/3+1 



Tj ~ C 

— — — r(/3 + 1) — — 

Ma? 2 



+ 2-C 



(48) 



Once again, here we are looking for polynomial solutions for the function G(x). In order to do that 
we must impose conditions equivalent to (1331) and (1341) on the coefficients. 
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From the condition C - 2n, we obtain the energy spectrum: 



k 2 + 



ni 



MV 2m w „ w / |m| , 
Mw + 2Mw « + — + 1 

f + M 2 w 2 a \ a 



(49) 



with u) - jj yjr]^ + M 2 io 2 . Like in the previous analysis, we can see that the energy spectrum of the 
particle depends on a given set of parameters. The condition d n +\ = provides an algebraic equation 
involving this set of parameters and the mass of the particle. Here, in this subsection, we shall adopt 
a different procedure as we did in the last subsection. We shall admit that the strength of the magnetic 
field is adjustable according this equation which in its turn will depend on the quantum numbers n and 
m. This will provide discrete values for Bo, and consequently for co which we define now by w n>m . 

To illustrate, let us apply this formalism for the case where the function G(x) is a polynomial of first 
order. In this case the magnitude of the magnetic field will be given by u)\^ m as exhibited in d52l ). The 
function G\{r) reads, 



G hm (r) = 1 + -^-r . 



(50) 



As to the self-energy we have: 



Elcjnl — i 



k 2 + 



M 4 a> 2 
T] 2 L + M 2 U) 2 



l.m 



2m _ ( \m\ 

— Moj\ m + 2Mu>\ m I — +2 
a ' \ a 



(51) 



with the parameter a)\ >m assuming the form 



rjt(\m\ + 3 
M \ a + 2 



(52) 



For any values of the physical quantities, equation (1521 presents at least one real physical solution for Bq 



4.. 4 General Solution 

In this subsection we shall investigate the solution of the problem under consideration in its more general 
form, which is the analysis of the quantum motion of a bosonic charged particle in the presence of a 
uniform magnetic field, given in terms of a vector potential by (O and a scalar potential given by © in 
the cosmic string spacetime. The corresponding radial differential equation associated with this problem 
is CD. 

As we did previously, in order to obtain the energy spectrum, it is convenient to analyze the 
asymptotic behavior of (fTT~|) for x — > and x — > oo. After we have made this analysis we can express 
the function u(x) as follows: 



u(x) = xP' 



e 1 



>H{x) , 



(53) 



with 



m- 



(54) 



Substituting the above function into (fTTT) . we get, 



xH" + 



2B+^x- 2x 2 



H' + 



H = 0. 



(55) 
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Once more, by using the Frobenius method, 



H{x) = Y t c n x n , 



(56) 



we can observe that the recurrence relations obeyed by the coefficients of this expansion is: 



1 



Cn+2 - 



(n + 2)(2/3 + n + 1) 



(- S ±(n+Y)-By n+l + (2n-C)c n 



1 (S L S c \ 



2 \ A j8 



(57) 



where 5 = /3^ + 6c, C = ^ -2/8- 1 + As in the previous cases, below we exhibit the first coefficients, 
considering cq = 1 : 



C2 



C3 



2y8+ 1 



/3Ma> 



6(J3+ 1)VMS Uy6+ 1 



4?7l 
— + 



+ M(2-C)H^ , (58) 



with O = ^ + ^ . To obtain a polynomial form of the function, we must impose two conditions on the 
coefficients similar to (1331 ) and (l34l ). The condition C = 2n, provides the energy spectrum, 



, MV m / 1 

r + — — + 2r] C 7] L - 2Mco— + 2Moj « + -+/? 



rj 2 L + M 2 u) 2 



a 



(59) 



being co = ^ + M 2 co 2 . We can see that (l59l does not change by the discrete symmetry ?/c — > -??c 
together with tj l — > -rj L . 

From this result we can see that the magnetic field and both terms in the scalar potential interfere 
significantly in the energy spectrum. For the case were the function H(x) becomes a Heun polynomial 
of first degree, 



the energy of the particle is 

Ek,m,\ ~ — 



k 2 + 



H hm (r) = 1 + M 



Ma) 



1 ,m 



P 



r , 



(60) 



2m 
a 



, -, +2t] C Vl Moj hm + 2Mco hm \fi + - 



(61) 



The condition C2 - imposes that the parameter oj\^ m satisfies the following third degree algebraic 
equation, 



Mr/, 



CO 



C,-,2 



l,m 



CO 



l,m 



■coi. 



M 



= 0. 



(62) 



Although the above equation presents at least one real solution, we decide do not reproduce it here 
because the latter is a very long expression. We also would like to emphasize that the energy level (loTT ). 
has been obtained assuming that the parameters rjc and rji are different from zero, C = 2 and C2 = 0. This 
is the lowest energy level that the analytical method adopted can provide. We can see that this energy 
level is completely different from the corresponding one obtained from the Ladau levels given by d23T ). 
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5. Fermionic Fields: Dirac Equation 



In this section we shall consider the relativistic quantum analysis of a spin-^ charged particle considering 
the most general situation, i.e., the fermionic particle in the presence of a uniform magnetic field and 
scalar potential, given, respectively by © and (f3]>. 

Let us first introduce the generalization of the Dirac equation to an arbitrary curved spacetime, in the 
presence of an external electromagnetic potential and a scalar potential, which is given by 

\i-f(x) (V„ + ieA„) ~(M + S (r))] - . (63) 

Here M is the mass of the particle and the covariant derivative for fermion fields is defined as follows 



together with the spinor afhne connection 



(64) 



(65) 



where is the Christoffel symbol of the second kind and e^(x) the tetrad basis defined later. Also, 
y /i (x) represents the generalized Dirac matrix satisfying the Clifford algebra 



{y"(x),y v (x)}=2^(x), 
and defined in terms of a set of tetrad fields ^(x) and constant Dirac matrices, y^ a) , as 

A*) = «f fl )(jf)7(a) ■ 
The tetrad fields, eft Ax), satisfy the relation 



rj^ b >^ a) (x)e^x)=g 



(66) 



(67) 



(68) 



Here Greek letters are used for tensor indices and Latin letters (in parenthesis) for tetrad indices. The 
matrices y^ are the standard flat spacetime Dirac matrices. In this paper we shall use the definition of 
such matrices, as given below: 



y 



(0) 



1 
-1 



and y 



(a) 







cr 



Ui) 







(69) 



where cr^ a) are Pauli matrices and 1 the identity 2x2 matrix. 

In order to consider the Dirac equation in the spacetime of a cosmic string, let us consider the tetrad 

e?, . as: 

(a) 



(a) 



1 











COS (f> 


_sin0 q 
ar 





sin0 


COS (f> q 

ar 


. 





1 . 


we 


obtain 


the followin 



(70) 



matrices y^(x): 



y° = y (0) , 



and y ! 



a 1 

-<T l 



(71) 
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Here, represents the following set of modified Pauli matrices: 

e~ H 



o~ — O" 



e i<P 
i 

ar 



cr 3 =cr = cr^\ 





e 



-er 







(72) 

(73) 
(74) 



hence y 3 = y^ 3) . 



5..1 Fermionic Radial Equations 

In this subsection, we shall investigate Dirac's equation d63l >. decoupling the differential equations for 
the four components of the Dirac spinor, through a series of symmetry considerations. As in previous 
sections, the presence of the magnetic field and the scalar potential is considered, just as described in 
section |2~1 

From equations (l65l ). (1701 and (1711 ). it is possible to evaluate the spinor connection 



r M = (o,o f i> > o), with r> = ^-^ 3 ), 



and T 



(3) 



O®* 




The Hamiltonian associated with (1631) reads 

;(i~< 



*H = -ij(o) 



y r d r + y nds + 



I (3) + ieA^r)] + y 0) d z + i(M + S(r)) 



(75) 



(76) 



This operator satisfies the commutation relations £3] = 73] = 0, where = -id z , J3 = L3+S3 = 
-id,/, + 5^ (3) , and the following linear eigenvalue equations hold 



= kW , 

J^V = -idfl + ^X (3), F = fP , 

where j = m + A, m = 0, ±1, ±2, . . ., and k € [-00, 00]. 

We shall adopt the following ansatz for the wave function: 



(77) 
(78) 

(79) 



-\Jr 



9+(r) _ 
-i<P-(r)e" 

X+(r) 
iX-(f)eP 



(80) 



and the radial functions obey the following coupled differential equations: 

/- + (—- Mcor) X - +k X+ + [M-E + S(r)] cp + =0, 
\ar I 

y - (A - Majr) X+ +k X - + [M-E + S(r)] <p-=0, 

<p'_ + (—- Mcor)(p- -ktp + + [M + E + S(r)]x+ = , 
\ar I 

<p'-(—- Mcur)<p + -ktp^ + [M + E + S(r)] y_ = . 
\ar I 



(81) 
(82) 
(83) 
(84) 
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Those equations exhibit the discrete symmetry 



X+ - Atp+ and if- = X\- > with A = - and s = ± 1 



(85) 



We have also introduced the cyclotron frequency as a> = jj^. Applying the symmetry argument expressed 
in (f85l) . we are able to reduce the number of equations to the only two given below: 



d 2 


j(j ~ <*) 


dr 2 


a 2 r 2 


'd 2 


Kj + a) 


dr 2 


a 2 r 2 



+ E 2 -k 2 + 2Mco 



+ E 2 -k 2 + 2Mco 



H 1 " 



M 2 oj 2 r 2 ~(M + S(r)f 



M 2 co 2 r 2 - (M + S(r)Y 



* + ^* = 0, (86) 
dr 

* + ^ = 0, (87) 
dr 



where we have used <p + = <p, and \- = X f° r simplicity 

Equations (l86l ) and (I87T ) are coupled second order ordinary differential equations. We notice that this 
pair of equations can be written in the following way 



K 



0, 



where K is the 2x2 matrix operator below, 

K = D(r)l + A(r)cr (3) + S'{r)cr (l) 



with 



d 2 i 2 i 
D(r) = — - -!— +E 2 -k 2 + 2Mco- - M 2 io 2 r 2 - (M + S(r)) 2 
dr 2 a L r l a 

A(r) = Moj , 

ar 2 



S (r) = — T + t]l ■ 



(88) 

(89) 

(90) 

(91) 
(92) 



The operator K can be diagonalized by using a 2 x 2 orthogonal matrix R as shown, 

0, with R 



a. a >v 1 -A- I LO 

RKR~ R\ * 



X 



cos(0/2) sin(0/2) 
- sin(6>/2) cos(0/2) 



This process allows us to write the new operator, named Kg, by 



Kg = D(r)l + (-^-+Mw 

ar 2 



cos y - sin f 
- sin 9 - cos ( 



sin y cos y 
cos - sin ( 



(93) 



(94) 



In conclusion, this differential operator is diagonal if both conditions below are simultaneously fulfilled: 



— tan 9 = —rjc and Moj tan Q - t]l ■ 

a 



This analysis produces the following set of relations 



to = - 



sin 9 = - 



cos 9 



m 

aMrjc 
aye 

yjj 2 + a 2 i] 2 c 



^j 2 + a 2 T] 2 c 



(95) 

(96) 
(97) 

(98) 
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This means that in order to decouple ((861 ) and dBTT ), one must impose the condition d96l ). Accepting this 
condition the analysis of the quantum motion becomes enormously simplified. The two fourth order 
radial differential equations obtained in the general situation, reduce themselves to two second order 
differential equations. Although this condition provides a much simpler system, in principle the solutions 
for the wave equations and energy spectra may not be compatible. Fortunately, as we shall see later, this 
situation does not occur in this analysis. Now returning to d96l ), we see that the magnetic field strength, 
the mass of the particle and the scalar coupling constants must be related, and, on the other hand, this 
relations depends on the total angular momentum quantum number, j. In order to fulfill this condition, 
we may choose that the parameters a> or tjl, must be adjusted to a specific value if we are inclined to find 
exact solutions to the radial equations. As we did in the bosonic case, and accepting the condition (|96l >, 
it is possible to rewrite the previous differential equations in a Schrodinger form, as follows: 



d 2 if 
dr 2 

fx 

dr 2 



v l A£ 2 i2^ n ,, j n til 7(7 - 1) 2Mt] C 2 7 ?l 2 

E - M - k + 2Ma> 2t]cT]l ~ 7 ~ 2Mr\ir - y —r 

a T] C r 1 r rf 



v2 ,,2 j,.,, Jr. , riL 7(7+1) 2Mt] C 
E - M — k + 2Moj 2r\cT]L + 7 5 2Mrj L r - y *—~r 



a 



Vc 



2 rfi ,,2 

4 



<p = 0, (99) 

x = o, 

(100) 



where 7 = ^ -^j 2 + a 2 r] 2 c . Moreover, we can write these equations as BCH ones given by (1 1 27b - 
Adopting rjc as a positive constant, and defining a new dimensionless variable 




(101) 



we have 



(p"(x) + 



ZkjVc tjl 7(7-1) a 2 
— — - — — + - + bx - x 

ym\ m\ * 2 x 

£kjnc r] L y(y + 1) a 2 

1 1 v bx - x 

7\t]l\ \t]l\ x 1 x 



<p(x) - , 
Xix) = , 



(102) 
(103) 



where 



ft i - E 2 - M 2 - k 2 + 2M(x>- - 2t] C t] L , 

a 



a = -2Mr] C 



arjL 



7\t]l\ ' 



7\Vl\ 



(104) 
(105) 
(106) 



We stop at this point and leave the obtainment of analytical solutions to (11021) and (11031 ) for section |6~1 



6. Analytical Solutions for the Dirac Equation of Motion 

Following in the same direction as we did to obtain solutions of the Klein-Gordon equation, it is 
convenient to analyze the asymptotic behavior of the solutions of the radial differential equations (1991 
and d 100b for r — > and r — > 00. Through this analysis, we verify that the radial components of the 
Dirac spinor can be expressed as 

cp(x) = x y e-^ x - b) F{x) , (107) 
X(x) = x r+1 e-^ x - b) H(x) . (108) 
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Inserting this ansatz into 01O2I ) and ( 11031 ), we obtain 



xF'\x) + {ly + bx- 2x 2 ) F\x) + |l + ^ J a 



+ 



—TT + ~T ^ -* _ i — ' x 



F(x) = , 



(109) 



xH"(x) + (2(y + 1) + bx - 2x 2 ) H'(x) + ( 1 + (y+ 1)??L ) 



+ 



~T~T + — - 2(y + 1) - 1 + — x 



tI?7lI 
#(x) - . 



These are BCH equations with solutions 



TO = H B \2y- 1, 6, — — - + — - : — r 

jVil\ 4 \?] L \ 

= 2y+l, Z>, ^^ + -r + ^ 
rl^Ll 4 |?7 L | 



2a; -x 
2a; -x 



(HO) 

(HI) 
(112) 



It is useful to follow the analyses of both equations d 1091) and dllOl ) simultaneously. To accomplish 
this we may define the constants A, 5 and D, and the function G representing F or H as shown below: 



A - 



7, if G - F 
7+1, if G = # 

and D = 



1 + g, if G = F 



1 + 



(y+i>7L 
ytozJ 



, if G = H 



^ + ^-2(7+1) + !-^, if G = F 
|f + ^-2(7+1)-! + ^, if G - H 



(113) 
(114) 



By making use of theses definitions, it is possible to write the two differential equations in just one single 
form, 



xG" + (2A + bx - 2x 2 )g' + {Ba + Dx)G = . 
Here we shall use Frobenius method to find series expansions to this solution: 

G(x) = J] C n x n . 

n 

Substituting the above series into d 1 151 ), we obtain the recurrence relation 

-1 



(115) 



(116) 



C n+ 2 



(n + 2)(n + 2A + 1) 



{[Ba + bin + 1)] C„+i + (D - 2n) C n \ 



Ba 

and Ci = -^Co ■ 



(117) 
(118) 



With this relation we can construct the first three coefficients of the expansion. Assuming Co = 1 they 
are: 



C 2 = 
C 3 = 



1 



Ba 

— (Ba + b)-D 
2A 



2(2A + 1) 

1 ( Ba + 2b 
6(A + 1) \2(2A + 1) 



D- — {Ba + b) 

2A V ' 



(119) 
(120) 



At this point we know that by breaking the series expansion of the BCH function into a Heun polynomial 
of degree n, we obtain an analytical solution to the radial equations. This can be implemented by the 
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imposition of two conditions on the coefficients in a way similar to what we did in the bosonic analyses. 
These conditions are: 



C„+i = and D - 2n , with n — \,2, ... 



(121) 



From the condition D = 2n, it is possible to obtain a formal expression to the energy. We should be able 
to reproduce the same energy eigenvalue independently of the choice we make for the definition of the 
constant D (either for G = F and G = H). For this to occur, we need to impose rji > 0, otherwise, there 
are no eigenvalues associated with the Heun polynomials which satisfy the condition D = 2n for both 
components. After adopting this restriction, we obtain: 



E kdn = ±\M 2 -k 2 - 2Moj- + 27] C j] Ln + ?] Ln 

a 



2y(y + 



M 2 7] 2 C 
7 2 Vl„ 



(122) 



If we assume the solution as a Heun polynomial of degree n = 1 , the equation C2 = provides that the 
parameter t]l will depend on the total angular momentum quantum number j. So for this case we have 



2M 2 r] 3 c (2y + 1) 



(123) 



The above value together with (11221) provide the energy spectrum of the particle. 



7. Concluding Remarks 

In this paper we have investigated the relativistic quantum motion of charged spin-0 and spin- \ particles 
in a cosmic string spacetime, in the presence of a uniform magnetic field parallel to the string and scalar 
potentials whose centers are on the string. The complete analysis is twofold: it is applied to bosonic 
particles in the first part followed by the fermionic case. 

For the bosonic case we obtained from the Klein-Gordon equation, a radial differential equation (|TT]) 
by adopting the ansatz © for the wave-function. The former is analyzed considering different values 
for the physical parameters. Specifically for the most general case investigated in subsection I4..4L we 
were able to express the corresponding differential equation in terms of the biconfiuent Heun differential 
equation. Because we are interested in obtaining an explicit expression for the energy spectrum, we 
searched for polynomial solutions for this radial equation. So, in order to obtain these polynomials, we 
used the Frobenius method and a three-term recurrence relation for the coefficients of the expansion was 
obtained. Finally, by imposing specific conditions on these coefficients, polynomials solutions are found. 
Admitting a first order polynomial solution, the energy found was given by (l6Tb under the condition that 
the cyclotron frequency is a real solution of (l62l . Because this solution is a long expression, we decided 
to not reproduce it in the paper. 

For the fermionic case, we observed that following the particular choice of the ansatz for the wave- 
function (l80l ). a set of four radial differential equations is obtained from the Dirac equation. Moreover, 
by adopting the discrete symmetry given by (I85T ). we were able to reduce the system to only two coupled 
differential equations. A direct attempt to decouple these equations would result in a much more complex 
system of fourth-order differential equations. Instead, through the relations (|951>-(|98T>. which require that 
both scalar potentials be considered (this is particularly clear in d95l)), it is possible to greatly simplify 
this set of differential equations into biconfiuent Heun equations. As happened in the bosonic case, it is 
also useful to simplify the BCH solutions into Heun polynomials in order to obtain the energy spectrum 
of the quantum particle. This provides an additional algebraic relation between the parameters of the 
system. 

In conclusion, we notice that a simpler analysis of the motion of the fermionic field can be obtained, 
provided that two conditions involving the parameters of the system be satisfied, namely (l96l > and (11231 ) 
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for the case where the Heun polynomial is of the first degree. Because we have accepted r\ c > 0, to 
obtain a physically relevant result for the energy spectrum, we have also to assume that > 0. As a 
consequence, in order to (|96l > be satisfied, for a magnetic field along the positive direction of z-axis, i.e., 
B - Bq1c, the value of the total angular momentum quantum number, j, must be negative. Having chosen 
the opposite direction to the magnetic field only positive values for j would be accepted. We conclude 
that the general case must be dealt with by means of two different analyses: for a magnetic field along the 
direction k a negative value of the total angular quantum number is required; moreover, for the magnetic 
field in the opposite direction, a positive value is necessary. The fact that the sign of the quantum number 
j is conditioned to the direction of the magnetic field does not represent a strong restriction in our result, 
since the energy spectrum obtained in (11221 ) depends on the product between u and j which is always 
negative. 

As we have seen, closed analyses involving Heun biconfluent equations can only be carried out by 
imposing a specific Heun polynomial of degree n for the solutions of the field equations. The reason for 
this difficulty is that, for the complete biconfluent Heun function, there is no general expression for its 
asymptotic behavior for large values of the argument]! In this present analysis, for bosonic and fermionic 
particles, we have illustrated how relevant physical informations about the systems can be obtained when 
the BCH function is expressed as a Heun polynomial of degree n — 1 (of course this analysis can be 
extended for any polynomial of degree greater than unity). As consequence of this approach, some 
physical parameter must be fixed by the order of the polynomial, n, and also by the angular quantum 
numbers. The analysis of a BCH differential equation was also developed by Vercin [26] and Vercin at 
al En . investigating the quantum planar motion of two anions with Coulomb interaction between them 
and in the presence of an external uniform magnetic field. In the first paper the author found exactly the 
energy spectrum, by imposing that the magnitude of the magnetic field obeys a closed relation with the 
angular quantum number. In the second, the energy spectrum was evaluated numerically. 

The Coulomb-type scalar potential presented in this paper, can be associated with the presence of an 
electrostatic self-interaction on the charged particle induced by the non-trivial topology of the conical 
spacetime. As it was shown numerically in |9), this self-interaction is positive for a = 1 - 4fi < 1, 
which corresponds positive linear mass density of the string //; an attractive self-interaction takes place 
for a negative value of fi. As to the linear scalar potential, it can be associated with an extra cylindrical 
harmonic oscillator acting on the particle, whose classical angular frequency, a) c , being given by tjl/M. 
Moreover, the linear potential, whose original motivation comes from the confining model for quarks 
adopted in (TTJ ffifl . was basically included in this paper to make the present analysis as general as 
possible. 

Due to the magnetic field and the linear scalar potential, only bound states solutions can be presented 
in the systems considered. By choosing specific relations between the scalar potentials strength and 
the magnetic filed, the analysis of the energy-spectra become drastically simplified. Accepting these 
relations we were able to obtain explicitly the spectra associated with both, bosonic and fermionic 
charged particles; moreover, from these results we can also observe that, although the quantum particles 
move in a region of spacetime that is locally flat, the correspond energy spectra of the particles depend 
on the global properties of the spacetime through the parameter a. In fact these energy levels depend 
on the the orbital angular quantum number for bosonic particles and total angular quantum number for 
fermionic ones in the combination m/a and j/a, respectively. These facts are consequence of the planar 
angle deficit associated with the conical structure of the two surface orthogonal to the cosmic string. 

Before to finish this paper we would like to mention two potential applications of our computation. 
They are in Condensed Matter and Astrophysics. As to condensed matter physics, it is well known 
that linear defect in elastic solid named disclination can be dealt with the same geometric approach as 
cosmic strings [28]. In fact the analysis developed here constitutes a relativistic extension, with some new 
ingredients, of the previous one presented in iTTTTl to study the quantum mechanical motion of charged 
particles under the influence of an external uniform magnetic field in an elastic medium with disclination. 

5 For specific values of the parameters, the BCH function can be expressed as the confluent hypergeometric function. 
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Moreover, specifically to the analysis in section [fTJ of the quantum mechanic motion of a fermionic 
charged particle, it can be adapted to study graphitic cones under the influence of uniform magnetic field 
within the framework of long-wavelength Dirac-like model for electronic states in graphene (for a review 
see IT29TD . With respect to astrophysical application, this analysis may also be useful to understand the 
quantum motion of charged particle in the cosmos under influence of galactic magnetic field considering 
the presence of a cosmic string^. Summarizing, we can say that, by the results found in this paper, the 
presence of a cosmic string produces significant modifications on the energy spectra associated with 
bosonic and fermionic charged particles the presence of uniform magnetic fields. 



A Heun's Differential Equation 



The Heun's equation is the generic differential equation with four regular singular points at 0, I, a and 
oo. The standard natural form of Heun's equation reads ETTl . 



y"k) + (- + -^-r + -^-W) 4 ,"?-' ; )*) - - 
\z z - 1 z-al \z(z - l)(z - a) j 



(124) 



The parameter e is expressed in terms of the other ones by e = a + cr+l- A-/u. Four confluent equations 
arise from the general Heun equation by means of different confluence processes. Of particular interest 
to our investigation is the biconfluent equation (BCH), with a regular singular point at z = and irregular 
singular point at z = °°. The canonical form of the BCH is 



xu"(x) + (1 + a - crx - 2x 2 )u'(x) + |(/1 - a - 2) x - -[^i + cr (I + a)]\ m \ ) = . 
where a, cr, A and fi are arbitrary parameters. By means of the transformation 

i+a trx+x 2 

u{x) = x 2 e 2 v{x) , 
the BCH equation takes on the Schrodinger form 



v"(x) + 



, D E 
Ax + Bx + C + h — 

x x 1 



v(x) - , 



(125) 



(126) 



(127) 



with A = -1; B = -cr; C = A- 2f ; D = -f; E = ±f-. 

The biconfluent Heun polynomial, Hs(a, cr, A,/u; x), can be obtained as a Frobenius solution to BCH 
equation computed as a power series expansion around the origin, 



(128) 



n=0 



Substituting the above expansion into (1 125b . a three-term recurrence relation for the coefficients is 
obtained: 



1 



im+2 



for m > 1 and 



(n + 2)(n + 2 + a) 



- + cr In + 1 H 



C„ + i + (2n - A + a + 2) C„ \ , 



A, = - 



cr + 



1 + a 



(129) 



(130) 



6 An estimated strength for our galactic magnetic field can be found in 1 30 1 . This value is of order 3 x 10 6 gauss. 



18 



The roots of A n+ \ are the eigenvalues corresponding to these particular Heun polynomials. In this case, 
all subsequent coefficients cancel and the series results in a polynomial form of degree n for Hg- For this 
case we must have, 



A-a-2 = 2n, for n = 0,1,2,3... and A„ +1 = . (131) 

Unfortunately, there is no closed expression for the asymptotic behavior of the biconfluent Heun 
function in the general case for large values of its argument. Consequently, the conditions imposed on 
its parameters in order to provide a polynomial form for them, can only be obtained by analyzing each 
coefficient separately as we did along this paper. 

The biconfluent Heun function can be expressed in terms of the Whittaker M function, or confluent 
hypergeometric function, for specific values of the parameters, as shown below: 

1 x 2 (A a 9 \ (a A 1 a 9 \ 

H B (a,0,A,0;x) = -^e-M[-,-;x 2] J = i F l + - - + 1; x 2 \ . (132) 
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